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MODELING VIRAL EPIDEMIOLOGY IN CONNECTED NETWORKS* 

WILLIAM M. SPEARS 1 ’*, LORA BILLINGS 2 ’*, AND IRA B. SCHWARTZ 2 
l Al CENTER - CODE 5515, 

Special project in nonlinear science - code 6700.3 

PLASMA PHYSICS DIVISION 

NAVAL RESEARCH LABORATORY, WASHINGTON, DC 20375 
February 14, 2001 

Abstract. We derive some Markovian and differential equation models of viral epidemiology 
on connected networks. We examine the stability properties of endemic states for models based on 
contact rates derived from probabilistic concepts of connectivity. In particular, boundaries which 
delineate the onset of viral propagation in a connected network are derived and analyzed. 


Key words, viral epidemiology, Markov chains, transient behavior. 


1. Introduction. We are becoming increasingly dependent on large intercon¬ 
nected networks for the control of our resources, such as the Internet, communications 
networks, and power grids. The advantage of these networks is the ability to route 
resources in a reasonably optimal fashion. However, their interconnectivity, coupled 
with the lack of global view of what is happening in these networks, can lead to 
tremendous problems in network reliability. For example, small local failures can eas¬ 
ily propagate to entire networks, causing loss of service and corruption of data. Also, 
deliberate attacks, such as viral attacks, can easily cause widespread havoc. The pur¬ 
pose of this paper is to outline a broad framework for describing the spread of viruses. 
The framework is designed to be as general as possible, with restrictive assumptions 
only introduced when needed. The focus will be on two complementary analysis tech¬ 
niques - discrete-time Markov chain models and continuous-time differential equation 
models. We give a brief overview of Markov chains first. 

2. Overview of IVtarkov Chains. A discrete-time Markov chain is a dynamical 
system composed of S discrete states. At each time step, the Markov chain can change 
states. Let Xt be the random variable for the Markov chain, which can take on any 
of the <S states at time t. The system is described by an <S x S matrix Q , which gives 
the probability of transitioning from one state a at time t — 1 to state b at time £, 
where a component of Q is given by: 


Q(a,b) = Pa, b - P{Xt = b | X t -i - a) 

The p a i values define the “one-step” probability transition matrix Q, since Q 
describes the probability of transitioning from state to state in one time step. The 
transient behavior of the system is obtained from the “n-step” probability transition 
values, which are obtained from the nth power of Q : 

Q n (a,6) = pa,b (n) = P(X t — b | X t ~ n = cl) 

It is also possible to compute conditional probabilities over a set of states. Define a 
predicate Precis and the set B of states that make Precis true. Then the probability 
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that the system will be in one of the states of B at time t , given it is in state a at 
time t — n is: 


P«,B (n) = P(X t £ B I X t - n = a) = J2 P«* (n) 

b€B 

Analysis of Markov chain behavior can be considered to be “instantaneous” or 
“cumulative.” Instantaneous behavior refers to events that occur at a particular 
time. An example would be whether the system is at a particular state at that time. 
Cumulative behavior refers to events that have occurred by a certain time. 

2.1. Instantaneous Transient Behavior at Time n. The probability of being 
in some state b at time n is given by simply considering the probability of each possible 
n-step transition, appropriately weighted by the a priori probabilities: 

P6 (n) = P(Xn = b) = 

a 

where the a priori probability of a system being in state a at time 0 is denoted as 
Pa (0) - 

The probability that the system is in one of the states of B at time n is: 

PB (n) = E p b W (2.1) 

Finally, it is possible to compute the probability that the system will transition 
from one set of states to another set of states. Let Pred a be another predicate over 
the states, and denote A to be the set of states that make Pred,A true. Then the 
probability that the system will be in one of the states of B at time t, given that it is 
in one of the states of A at time t - n, is: 

PA , B {n) = (2.2) 

P(X t G B | A f _n G A) = n> Pa - g(n) _ L) a€ ytPa (< n) Pa,B (n) 

p^ (t - n) £. € a p« (t - n) 

which involves a renormalization over the states indexed by A. Since Eq. 2.2 describes 
how a system transitions from a group of states to another group of states, this 
raises the intriguing notion that a system with a large number of states might be 
simplified (or aggregated) into a system with a smaller number of groups of states. 
Unfortunately, this is hard to do in general since Eq. 2.2 is a time-dependent equation 
(thus making the aggregated Markov chain non-stationary), however some progress 
has been made in this area [3]. The nice feature of this formalization is that any 
predicate over the states can be used. 

It is possible to generalize further to arbitrary functions / over the states and 
compute, for example, the expected value of that function, at time n: 

E[f} (n) = EP^ (n) /W (2.3) 

b 

2.2. Cumulative Transient Behavior. Another common computation involv¬ 
ing Markov chains is referred to as the “mean first passage times” for going from state 
a to state b (for a nice discussion of this, see [4]). This refers to the length of time that 
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it takes (on the average) to reach state 6 for the first time, given that the process is 
currently in state a. Answering such questions involves solving the set of simultaneous 
equations: 


m a ,b = Pa,6 + Y JV* ( 1 + m M) ( 2 - 4 ) 

k^b 

where m a ,b denotes the mean first passage time from state a to state b . To understand 
the equation, consider transitioning from state a to 6 in one move. This occurs with 
probability p a ^ and requires only one step. However, suppose the system transitions 
from state a to state /c, where k is not equal to b. This occurs with probability p a ^ 
and requires one step. However, there now remain steps to state b. 

As before, if there is interest in a set B of states, it is possible to compute the 
mean first passage time for the system to first enter that set of states, given that it is 
currently outside that set: 

Pa,b 4- (1 +rafc,£?) 

b€B k£B 

where m a ,B denotes the mean first passage time from state a to any of the states in 
set B , and a is not in B. This is very similar to Eq. 2.4, with the exception that the 
probability of entering state B in one step is simply the sum of the probabilities of 
entering each state within B. 

Once this system of simultaneous equations is solved, it is possible to calculate 
the “expected waiting time” to reach a state in B , given a random initial state, via: 

EWTb = Y P« (0) 0 + Y Po (0) m a,B 

a£B a&B 

There are two parts to this equation. The first part reflects the possibility that 
a random initial population is in state B , and hence has a zero waiting time. The 
second part reflects the mean passage time from initial populations not in B , to a 
state in B. Clearly this simplifies to: 

EWT b = Y P* (0) m «' B 

agB 

This holds for any set of states J9, and thus it can be used to provide expected 
waiting times for a variety of events. 

3. Markov Models of Viral Spread. In this paper we are concerned with the 
spread of viruses throughout a network. We assume that the network consists of N 
nodes. Each node can be in one of four medical conditions: Susceptible (S), Exposed 
(E), Infected (I), and Cured (C). It is often customary to also refer to these conditions 
as “states”, but that would cause some confusion with the Markov definition of state. 
In this paper the number of nodes in the various medical conditions will define a 
Markov state. In general, if there are m medical conditions and N nodes, the number 
of Markov states S can be shown to be: 

•^("rr 1 ) 
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3.1. The Two-Condition I-C Case. We first start by ignoring the Susceptible 
and Exposed conditions, and assume that the nodes can either be infected or not 
infected (cured). Infected nodes can remain infected or become cured. Likewise, 
cured nodes can remain cured or become infected. This is the simplest case to analyze. 
Figure 3.1 illustrates the transitions between these two medical conditions, with their 
associated probabilities. 


5 



V' 

Fig. 3.1. The transition diagram showing how a node can change medical conditions. 


Let I be the number of infected nodes and C be the number of cured nodes. 
Since each node must be infected or cured at any given time, / + C = N. The 
Markovian state is defined by the pair (/, C) and we want to compute the probability 
of transitioning from state a = (/, C) to state b = (/',C'): 

Pa,b = P(X t = (I' 9 C) \X t - X =(/,C)) 

Such a transition means that V — I more cured nodes are infected than infected 
nodes are cured. Note that if V < / then the previous sentence is still mathematically 
and grammatically correct, but it is perhaps more understandable to say that I — V 
less cured nodes are infected than infected nodes are cured. 

Let x be the number of infected nodes that are cured. Then since there are I 
infected nodes, I — x of them are not cured. If we use 8 to denote the probability 
of curing an infected node, then the probability of curing x of the infected nodes is 
simply: 


mi-*) 


I~x 


Also, since V — I more cured nodes are infected than infected nodes are cured, 
we know that x *F V — I of the cured nodes must become infected. Since there are C 
cured nodes then C — (x + I* — I) of the cured nodes are not infected. If we use \j! 
to denote the probability of infecting a cured node, then the probability of infecting 
x + V — I of the cured nodes is simply: 

( x+ p _ j ) tt) x+, '-\ i - 

We can put this all together by summing over all possible values of x: 

p “- b = y 
X v ' 




c 

x + r -i 


(3.1) 


It now only remains to bound x correctly, which can be achieved by examining the 
combinatorials. The first combinatorial implies that 0 < x < L The second implies 
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that I - I* <x <C - V + /. In order for all of these constraints to be held we get: 


mtn{/,C— /' + /} 


Pa,6 


£ (0 


(*)*(! 


(3-2) 


(*+£-/ ) ( '‘ T+ '’"' (1 


Since C = N - I, this equation can be slightly simplified, but we do not do that 
here, in order to provide a smooth transition to the three-condition case. It is also 
useful to note that the number of Markov states for the two-medical condition case is 
simply: 


5 = 


N + 2-1 
2-1 


J = N + l 


3.2. The Three-Condition S-I-C Case. We now add in the Susceptible con¬ 
dition. Susceptible nodes can remain susceptible or become infected. Infected nodes 
can remain infected or become cured. Cured nodes can remain cured or become sus¬ 
ceptible. Figure 3.2 illustrates the transitions between these three medical conditions, 
with their associated probabilities. 
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Fig. 3.2. The transition diagram showing how a node can change medical conditions. 

Again, let / be the number of infected nodes and C be the number of cured nodes. 
Let S be the number of susceptible nodes. Since each node must be susceptible, 
infected or cured at any given time, S + / + C = N. The Markovian state is defined 
by the triple ( S , /, C) and we want to compute the probability of transitioning from 
state a = ( S,I,C ) to state b = ( S',I',C'): 

p a<b = P(x t = {S'J',C') I x t -i = ( S,r,c )) 

Such a transition means that I' - I more susceptible nodes are infected than in¬ 
fected nodes are cured. It also means that S'-S more cured nodes become susceptible 
than susceptible nodes become infected. 

Let x be the number of infected nodes that are cured. Then since there are / 
infected nodes, I — x of them are not cured. If we use 6 to denote the probability 
of curing an infected node, then the probability of curing x of the infected nodes is 
again: 
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Also, since V — I more susceptible nodes are infected than infected nodes are 
cured, we know that x + V — I of the susceptible nodes must become infected. Since 
there are S susceptible nodes then S — (x 4- /' — /) of the susceptible nodes are not 
infected. If we use /z' to denote the probability of infecting a susceptible node, then 
the probability of infecting x + V — I of the susceptible nodes is simply: 

( x + r- r ) (^T +r -'(i - / o s - (x+r “ ,) 

Finally, since S' -S more cured nodes become susceptible than susceptible nodes 
become infected, we know that x 1- I f — I + S' — S of the cured nodes must become 
susceptible. Since there are C cured nodes then C - (x +1' - I + S' - S) of the cured 
nodes are not made susceptible. If we use 5' to denote the probability of making a 
cured node susceptible, then the probability of making x + I' — I -F S' — S of the cured 
nodes susceptible is simply: 

f C ^ /r/\ z+r-i+s'-S(i _ r ,,c-( x +t-i+s'-s) 

V x + r-i + s 9 -s) { J [ } 

Putting this all together we get: 


mtn{/,C-J' + I-S'+S} 
Pa,b = Y, 

x=max{OJ—l* J — F'+S—S '} 




x + r-r) - A ,o 5_(x+r " ,) x 


(3.3) 


(s + r-? +S 


C-(x+I'-I+S'-S) 


For the three-medical condition case the number of Markov states is given by: 

0 ( N + 3-1 ^ _ (N + 2)(N + 1) 

S ~\ 3-1 )- 2 


3.3. The Four-Condition S-E-I-C Case. It is possible to add in the Exposed 
condition. A susceptible node becomes exposed before it becomes infected. Figure 3.3 
illustrates the transitions between these four medical conditions, with their associated 
probabilities. 



Fig. 3.3. The transition diagram showing how a node can change medical conditions. 
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The Markovian state is defined by the tuple (S, E, I,C) and we want to compute 
the probability of transitioning from state a = (S, E, /, C) to state b = (S',E ', C'): 

Pa, b = P(Xt = (S',E',I',C') I X t -i = (S,E,I,0) 

Such a transition means that I 1 -I more exposed nodes are infected than infected 
nodes are cured. It also means that E' - E more susceptible nodes become exposed 
than exposed nodes become infected. Finally, S' — S more cured nodes become 
susceptible than susceptible nodes become exposed. 

It is easy to produce p a , b by extending Equation 3.3: 


iin{I,E-I' + I,S-I' + I-E'+E,C-I'+I-E'+E-S'+S} . 

V—\ / i 


Pa, b = yz 

x=max{0,I-I'J-I'+E-E' .I-I'+E-E'+S-S') 


( I \ 
\ X J 


(sni-s) 


I-x . 


( , + r - f + f - «)-<-'■-—> x 

^ (( 5 , ) z + / '- / + £ '- £:+s '~ 5 (1 - tf} c -( x+r - r+B '~ B+S '~ s ) 


,l\x+l' — l I 


E-{*+I'~I ). 


C 


x + I' — I + E' - E + S' 


(3.4) 


For the four-medical condition case the number of Markov states is given by: 


3.4. Discussion. Given the above Markov models, we can use them to determine 
various quantities of interest. For example, the probability distribution of infected 
nodes at time n can be obtained, as well as the probability of extinction. The expected 
number of infected nodes at time n is easily computed. Finally, the mean passage 
time until viral extinction is easily obtained. 

Given these preliminary models, it should be clear at this point that they can 
easily be expanded to include m medical conditions, assuming that the transitions 
between medical conditions occur in a ring (as is true for the models examined thus 
far). Another important point is that the parameters 6 , p ', a, and 6' need not be 
constants - they in fact can be functions that depend on the state that the process 
is in. This level of generality can provide for nonlinear dynamics of the system, as 
we will see in the following discussion concerning differential equation models of viral 
spread. 

4. Differential Equation Models of Viral Spread. We can now provide 
analogous differential equation models. 

4.1. The Two-Condition I-C Case. If we let i = I/N and c = C/N be the 
proportion of infected and cured nodes (respectively) then: 


di , 

— — cu — 10 
dt ' 

dc , 

— = id - cu, 
dt 


(4.1) 
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Since i + c = 1 this can be simplified to: 
di 

— = c/i - iS = (1 - i)fi - id (4.2) 

where n' is the probability of infecting a cured node and 6 is the probability of curing 
an infected node. 

4.2. The Three-Condition S-I-C Case. If we add in the Susceptible condition 
then the system can be described as follows: 



where s = S/N is the proportion of susceptible nodes. In this case fj is the probability 
of infecting a susceptible node, J is the probability of curing an infected node, and 8 f 
is the probability of making a cured node susceptible. This could be reduced to two 
equations if desired, assuming a constraint of s -f i + c = 1. 

4.3. The Four-Condition S-E-I-C Case. If we add in the Exposed condition 
then the system can be described as follows: 


ds 

— — cd — sa 
at 


de 

; Tt =sa ~ e > J 


di 


dc 

dt 


i8 - c8‘ 


(4.4) 


where e = E/N is the proportion of exposed nodes. In this case // is the probability 
of infecting an exposed node, 8 is the probability of curing an infected node, 8 f is the 
probability of making a cured node susceptible, and a is the probability of making 
a susceptible node exposed. This could be reduced to three equations if desired, 
assuming a constraint ots-{-e + i + c = 1. 

4.4. Discussion. Again, it should now be very clear that the differential equa¬ 
tion models can easily be expanded to include m medical conditions, assuming that 
the transitions between medical conditions occur in a ring (as is true for the models 
examined thus far). Also, it is important to reiterate that the parameters 8 , fi *, a, 
and 8' need not be constants, but may be functions that depend on the state of the 
process (e.g., the proportion of infected nodes i). A concrete example is shown in the 
following section. 

5. Random Directed Graphs and Measures of Network Connectivity. 
The nice aspect of the above framework is that we have not had to yet define //, a, 
8 , and 8‘. In this paper we will assume that 8 is a constant that reflects an infected 
node’s ability to cure itself (the strength of the immune system or the ability of the 
sysop to rid the node of the virus). In the three-condition case, 8' is a constant 
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that reflects a cured node’s frequency of opening itself back up for infection (by 
weakening the immune system through malnutrition or becoming lazy as a sysop). 
In the four-condition case, a is a constant that reflects the probability of making a 
susceptible node exposed, and could be analogous to the likelihood of going to the 
mall or connecting a computer to the net. 1 On the other hand, the probability that 
a cured/susceptible/exposed node (depending on whether we are talking about the 
two-, three-, or four-condition case) will become infected, p 1 , is largely dependent 
on the number of other neighboring nodes in the graph that are also infected. The 
computation of p' is thus very dependent on the structure of the graph. In this paper 
we recalculate this quantity for the case of the random directed graphs defined in [1]. 

A random directed graph of N nodes is constructed by making random, indepen¬ 
dent decisions about whether to include each of the N(N - 1) possible directed edges 
that can connect two nodes. Each edge is included in the graph with probability p. 
Since there are N nodes, you can have at most N — 1 neighbors pointing towards each 
node (one doesn’t connect to oneself). In expectation any given node has p{N - 1) 
neighbors to itself. In [1] this is denoted as b = p(N - 1), and is a measure of the 
connectivity of the network. We now state the function for p' for three separate 
situations. 

5.1. Strong Links. In this case you can only become infected from your neigh- 
bors that are infected. The probability that an infected neighbor can transmit the 
infection to you is Thus the probability that an infected neighbor won’t transmit 
the infection to you is 1 —>9. Since i is the proportion of infected nodes in the network, 
in expectation you will have iNb/(N - 1) infected neighbors (because once again the 
assumption is that you yourself are not infected). The probability that none will 
infect you is (1 - Thus the probability that you will be infected by at 

least one infected neighbor is: 

// = i - (i - = i - (i - p) iNp 

where once again i is the proportion of infected nodes. 

If N is large, the reasonable approximation (which we will use throughout) is: 

// = !-(!- pf 


If y and z are sufficiently small it is well known that a reasonable approximation 
to —y) z is 1 — yz. If this approximation is taken (which is not necessarily reasonable 
in this situation) we obtain: 

// = Pib 

which agrees with [1]. Notice that although the systems of the two-, three-, and four- 
condition cases appear as linear differential equations, the effective contacts between 
infectives and susceptibles create a mass action term. This, in turn, makes the problem 
of viral spread in a network a nonlinear one. 

5.2. Weak Links. In this case you can also be infected from infected nodes that 
are not your immediate neighbor. The probabilit)' of transmission in this case is 
In expectation you will have N — b— I non-immediate neighbors, of which i(N -b- 1) 
are infected. These neighbors will not infect you with probability (1 - ^ U ,)* (W 6 l) . 


1 There is no reason why we can’t relax these assumptions if we wish. 
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Putting it all together, the probability that you will be infected by at least one of 
your immediate or non-immediate neighbors is: 

= pf (i-p u ,f N - b -" 

Making the same approximation used above, and throwing out further second- 
order terms (which again may not be reasonable to do) yields: 

//' = f3ib + /3 w i(N -6-1) 

which agrees with [1]. Very weak links in the form of mass action terms such as that 
derived above, can qualitatively change the dynamics in coupled populations [2]. 

5.3. Hierarchical Links. The final case we consider is that of hierarchical links. 
In this case you do not assume that you have 6 neighbors. Instead, neighbors are 
distributed around each node in L layers. The number of nodes at layer (distance) l 
is 2 l ~ l and there are a total of N = 2 L nodes. The probability of transmission from 
an infected node at distance l is given by j3 0 r*™ 1 , where r is a constant ranging from 
0 to 1. 

Following a similar process used above it is simple to show that the probability 
of being infected is given as: 

— 1 — II (! — A)7" i_1 ) i2l_I 
1 

Making the same approximation used above, and throwing out all higher-order 
terms (which again may not be reasonable to do) yields: 

L 

n' = J^A)(2 r)‘~ l i 

i=i 

which agrees with [1]. 

6. Analysis and Comparison of the Models. Now that /z' has been com¬ 
puted for random-directed graphs, we are directly able to use either the Markov chain 
or differential equation models of virus spread on random-directed graphs. We now 
run simulations of the Markov models and compare the results to the differential 
equation models. 

6.1. The Two-Condition I-C Case. The differential equation for the two- 
condition case is 

di — 

where the constant solutions are found by solving — = 0. The zero solution, i 0 — 0, 
exists for all f3 and 5. The other solution is called the endemic solution, which is 
implicitly defined on 0 < 8 < min(—61n(l —/?), 1). The two possible behaviors of this 
model are divided in parameter space by the curve 

S = -61n(l - P). 

If S > -61n(l - /?), solutions tend toward the solution z 0 for all initial conditions. In 
other words, the virus will naturally die out. In the rest of the parameter space, the 
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Fig. 6.1. The stability of io and i e in the two-condition case as a function of fi and 6. 



Fig. 6.2. The agreement of the Markov and ODE models in the two-condition case as a function 
of the number of nodes, N. In this computation, fi = 0.2, 6 - 0.2, 6 = OA(N-l), all of the nodes are 
initially infected, and the Markov model is iterated 20 time steps. The percentage of nodes infected 
increases exponentially to 0.83 percent. 


solution i 0 is unstable and all solutions tend toward the stable, endemic solution i e . 
Here, the virus will not naturally die out, but persist at a constant level. See Figure 
6.1. The preferred behavior is for io to be stable. 

The ODE model agrees with the Markov model depending on the initial condi¬ 
tions. Shown in Figure 6.2, as the number of nodes is increased, the percentage of 
nodes infected increases exponentially. The region where io is stable is also predicted 
by the Markov model, with the division matching the curve <5 = -61n(l - P) as shown 
in Figure 6.3. 

6.2. Tlie Three-Condition S-I-C Case. In two variables, assuming s+i + c = 
1, the differential equation for the three-condition case is 

dt (6.1) 

^U(l-(l -Pf b )s~6i. 
dt 

The zero fixed point is p 0 = (0, 1). The eigenvalues of the linearization about po are 
A, = -S' and A 2 = -5 - £ln(l - (i). Since 0 < S' < 1, then the eigenvalue A, < 0. 
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Fig. 6.3. For the two-condition case, the stable solution i predicted by the Markov model after 
20 time steps as a function of ft and 6. Here, N = 100 and b ~ 0A(N — 1). Superimposed in green 
is the curve 6 = —5ln(l — p) from the ODE model dividing where io and i e are stable . 


The other eigenvalue, A 2 , is negative for 1 > 6 > -61n(l - £), where p 0 a stable fixed 
point. Note that this is the same curve as in the two-condition case. The associated 
eigenvectors are 


ei =(0,1) 


e 2 


( 6-6' + bln(l-P) \ 
V <S'-61n(l-£) ’ ) 


( 6 . 2 ) 


In the rest of parameter space, p 0 is an unstable saddle fixed point. That is because 
0 < £ < min(—Mn(l — /?), 1), making the eigenvalue A 2 non-negative. 

The endemic fixed point, p e - ( i € ,s € ), is defined implicitly: 


0 - (1 - i € )S' - i e S - 

- 

1-(1 -0)*J 


i € SS ( 
1 - (1 - 


(6.3) 


For an example solution, see Figure 6.4. The eigenvalues for this fixed point are 
complex conjugates with norm less than 1, or stable, for parameters where the zero 
solution, po is unstable. The expression $ = —b ln(l—/?) defines a surface in parameter 
space (/?, 6 , (S') on which p 0 has neutral stability. This surface divides where the virus 
persists from where the virus dies out. 

It would now be possible to run the three-condition case simulation for the Markov 
model and compare the results to the ODE model, although we leave that for the 
future. 


7. Conclusions. We have derived a probabilistic framework that provides a 
more accurate representation of how viruses spread throughout a network. In par¬ 
ticular, the framework lays down the building blocks on which to incorporate more 
realistic parameters, such as measures of connectivity, effective contact rates, and 
different network architectures. 
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Fig. 6 .4. Sample solution in the three-condition case for initial condition (z,s,c) = (0.6,0.4,0), 
P = 0.1, 6 = 0.3, S' = 0.1, and 6 = 5. 

The framework is expressed in terms of Markov chain and differential equation 
models. Although the Markov chain models are inherently more accurate, they are 
also computationally more intractable. The differential equation models provide a vi¬ 
able alternative approach. Simple dynamic models expressing the state variables have 
been derived for large populations. The effective contact rates have been expressed in 
terms of network connectivity measures, and generate a viral spread that is the result 
of nonlinear mass action. 

Analysis of the models shows that for an endemic state to persist, the trivial 
solution of all susceptibles and no infectives must be unstable. A curve defining neutral 
stability between die out and endemicity has been derived and computed based on a 
stability exchange between the two states. Such curves provide a metric in which to 
measure the spread of a virus in a network, and may provide a future parameter in 
which to control the spread without sacrificing much network information flow. 
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